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ABSTRACT 

Time delays in gravitational lenses can be used to determine the Hubble constant and 
the lens potential. In future surveys, many gravitational lenses can be discovered, and 
their time delays and image positions can in principle be measured. Using an elliptical 
power-law potential, we show that combinations of image positions and time delays for 
quadruple lenses yield simple analytical expressions that are connected with observable 
quantities. These relations can be used to obtain the approximate axis ratio q, the 
Einstein radius and the slope. We apply this method to RX J1131—1231, and show 
that our analytical results match the full numerical determinations approximately. 
Our approach can quickly determine rough values of lens parameters, which can then 
be used as initial guesses for further refinement through numerical modelling and may 
be useful for automated lens search in large surveys. 

Key words: Gravitational leasing: strong - galaxies: structure - quasars: individual: 
RX J1131-1231 


1 INTRODUCTION 

Time delays in gravitational lenses provide an independent 
way of measuring the Hubble constant (Refsdal 1964). Many 
studies have been performed in this direction. Observation- 
ally, two dozens or so lenses have their time delays accurately 
measured after painstaking efforts (e.g., Tewes et al. 2013). 
Theoretical models also become increasingly sophisticated. 
For example, a recent study by Suyu et al. (2013) takes into 
account of the line of sight structures in addition to the lens- 
ing galaxy. Future surveys such as LSST will discover tens of 
thousands of lenses. In addition, hundreds of lensed super¬ 
novae will be found. A large fraction (~ 40%) will have crude 
time delays (Oguri & Marshall 2010), while a smaller frac¬ 
tion of these (~ 400) will have well-measured values (Liao 
et al. 2015). Such a large sample of lenses with time-delays 
can in principle be a very efficient way of constraining dark 
energy (Linder 2011). 

The time delay for an isothermal potential or density 
distribution does not depend on the angular profile of the 
lensing galaxy (Witt et al. 2000). Other than this, very few 
analytical results are known when the profile deviates from 
the isothermal shape. With thousands of time delays becom¬ 
ing available in the future, it is important to understand 
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better how the time delays depend on the density prohle, 
and their impact on the probability distribution of time de¬ 
lays. In this case one should not only focus on the time 
delays and the ratios of time delays measured in double and 
quadruple lenses. The ratio of time delays may depend on 
the lens model in a complex way. Therefore a diagram (of 
ratios) of time delays may be difficult to interpret. Many 
Monte Carlo simulations might be necessary to understand 
their behaviour and the dependency of different lens models 
on the parameters. In this paper we show that it can be use¬ 
ful to focus on the sum of time delays, and combinations of 
time delays and image positions. These quantities are more 
naturally connected to properties of the lensing galaxy such 
as the axis ratio or the slope of the inner potential. 

In this work we shall recall the properties of the ellipti¬ 
cal power law potential (§2). While this lens potential is not 
as realistic as the elliptical density distribution (e.g., Kor- 
mann et al. 1994; Keeton et al. 1997), it often resembles the 
latter when the ellipticity is small (Kassiola & Kovner 1993) 
and provides useful insights into the lensing properties. In 
§3, we introduce a new moment approach through which we 
derive new connections between the time delays and their 
image positions. We apply the results to the quadruple lens 
RX J1131 —1231 in section §4, followed by a short discussion 
in §5. For ease of reading, we put most of the mathematical 
details in the appendices. 
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2 LENSING PROPERTIES OF AN 

ELLIPTICAL POWER-LAW POTENTIAL 

The elliptical power-law potential widely used in gravita¬ 
tional lensing (e.g., Blandford & Kochanek 1987; Kochanek 
& Blandford 1987) is given by 

for 0 </3 < 2 (1) 

where {x, y) are the (angular) lens position, q is the axis ratio 
of the lensing galaxy (0 < g ^ 1), /? the slope of the potential 
and a sets the angular scale of the lens. For /3 = 1 we have the 
elliptical isothermal sphere which has very simple properties 
(see below). 


2.3 Time Delay 

The time delay in gravitational lenses is given in general by 
(e.g. Schneider et al. 1992) 

At = ^{1 + Zd)T{x,y), (8) 

2c 

T{x,y) = {x-0^ + {y-yf -2ij{x,y), (9) 

where Zd is the redshift of the lens, and D = DdDs/D^ b, 
with Dd,Ds and Dds are the angular-diameter distances to 
the lens, to the source, and from the lens to the source, 
respectively. Below we shall focus on r, and ignore all the 
cosmological dependences. 


2.1 Lens equation and Einstein radius 

The lens equation for the power-law potential is 

C = x-^=x-ax{x^ (2) 

^ ^ ( 3 ) 

where (^, y) are the source position. If the lens is circularly 
symmetric {q = 1) and perfectly aligned with the source, the 
images form an Einstein ring with its angular scale given by 

= (4) 

2.2 Image Positions, Numbers and Magnifications 


3 TIME DELAYS AND IMAGE POSITIONS IN 

QUADRUPLE LENSES 

As mentioned before, when a source is located inside the 
caustics, we may observe at least 4 images and possibly 3 
time delays if the lensed source varies. Such observations of¬ 
fer an exciting opportunity to probe the (inner) potential of 
the lensing galaxy. In particular, future surveys will observe 
thousands of gravitational lensed systems (see Oguri & Mar¬ 
shall 2010). Some of theses systems might be fairly isolated 
with negligible shear, which, as we show below, would be an 
ideal case to study the potential slope and the axis ratio of 
the lensing galaxy in more details. In any case, the effects 
of a (small) shear enter linearly (see §A1.3), and so can be 
ignored as a first approximation. 

For an isothermal potential with arbitrary angular pro¬ 
files, = rF{(j)), the time delay between two images 

depend only on the distances of the images {xi, yi) from the 
galaxy centre, namely 


The magnification of an image is generally given as 


1 

det J 


f d^dy 

d^dy\ ^ 

\dx dy 

dy dx J 


( 5 ) 


where {xi,yi) is the image position of the ith image in the 
lens plane. The image positions where the magnifications 
are infinite form critical curves and their source positions 
form caustics. For an elliptical potential, the caustics are 
in general diamond-shaped. When a source falls inside the 
diamond caustics, we can have either four images for 0 < 
/3 ^ 1 or five images for 1 < /3 < 2. In the latter case, a 
faint, central image appears. 

The total magnification is given as the sum of all abso¬ 
lute magnifications 






i 


( 6 ) 


For an elliptical power-law potential the sum of the signed 
magnification follows a simple relation 


2-/3 


( 7 ) 


for a quadruple system when the source is located inside the 
caustics. The result is exact for the cases /3 = 1 and /3 = 4/3 
(Witt & Mao 2000 and Hunter & Evans 2001). 


Ari,j = T{xi, yi) - T{xj,yj) = r| - r| (10) 

with Vi = \/xf + yf (cf. Witt et al. 2000). In particular for 
an isothermal elliptical potential it is rather straightforward 
to derive some invariants for the time delay. We can extend 
the results of Witt & Mao (2000) and obtain some simple 
relations concerning the moments of the magnification and 
image positions. Similar moment relations can be computed 
for the time delay as well. In particular one should not con¬ 
sider only the ratio of the time delays, but also the sum of 
all time delays of a system. We pursue this approach be¬ 
low; most of the mathematical details are collected in the 
Appendices. 


3.1 A Moment Approach 

We will develop here a moment approach using the time de¬ 
lays to study the deviation from the ideal case of an elliptical 
isothermal sphere. Let us denote the moments of the time 
delay by 

7 /"^ = "^TixT, ( 11 ) 

i 

Tym = 'Y^XiyT, ( 12 ) 

i 

where Ti = rixi, yi) is the time delay at the i-th image posi¬ 
tion. Since the absolute time delays n can not be measured 
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directly we relate them to measurable quantities by 


7/™ = ^Ari.ia;™ -f n xT 

(13) 

Tyn. = ^ Ari.lJ/r+T-l^yr 

i i 

(14) 

Now we use the Oth-moment equation 


To = ^ Ari,i -1- nn 

(15) 


i 


to eliminate the unknown time delay of the first image ri, 
where n is the total number of images, n is usually 4 (n = 5 
if a central image is present). Finally we can write 


AT.™ 

= 7- „ _ 2° a;’" 

n 

i 


and 

i i i 

(16) 

AT,™ 

i 



= ^Ari.iyr- i^Ari.i^yr 

(17) 


On the left side of the two previous equations we have now 
quantities which can be expressed for example in terms of 
77 , q and /3. On the right side, we have only terms involving 
the image positions and the relative time delays, which can 
be measured directly for a quadruple lens system. 


For fairly round system with g —>■ 1 the image positions 
are very close to the Einstein ring so that the solutions ap¬ 
proach ri —>■ 1 and the sum of the squared image positions 
^. rf —>■ 4. The time delay between the images becomes 
smaller and smaller so that the relative time delay moment 
vanishes (A 7).2 —>■ 0). For more elliptical systems the sum of 
the squared image distances to the galaxy centre increases 
with smaller q. The more elliptical the system becomes the 
more the image positions are located further away from the 
Einstein ring. In particular in clusters where one expects 
1 < /3 < 1.5 this effect is considerably amplified since for 
P = 4/3, B = 4 and for P = 3/2, B = Therefore in 
clusters the source location plays less of a role and the sys¬ 
tem is mostly determined by q and p. Similar effects can 
be observed if we have an isothermal sphere {p — 1) plus 
shear. The shear acts like an amplifier for smaller q. There¬ 
fore the shear 7 may emulate a more elliptical system than 
the ellipticity of the galaxy actually indicates. 

For the relative time delay moment A 7/2 the term be¬ 
comes small when the system is nearly round and q < 1. In 
this case the time delay is dominated by the source position 
and the approximation is fairly bad. For these cases we need 
to take the source positions into account. This would mean 
we need to solve a nonlinear set of equations and take the 
hrst order moments into account as well. Also in the case 
when the source is located near a cusp ^ or 7 might not 
be negligible any more. However, for rather elliptical sys¬ 
tem with q < 0.75 and p m 1 the source position plays less 
a role, so that the approximations hold fairly well in the 
case of galaxy leasing. Such a system would be suitable to 
determining the axis ratio q and the slope of the potential 


3.2 Analytical moment relations 


In the Appendix A we present the detailed procedure 
through which we find some analytic moment relations in¬ 
volving the relative time delay moments for a power-law po¬ 
tential. For simplicity, we normalise the coordinates by the 
Einstein radius, after which we have a — 1. For the second 
order moments we hnd, when the source is located close to 
the origin with no shear (7 = 0 ), 

2 +^ with (18) 


and 

A 7/2 = A%.2 + ATy2 ~ - ^ • (19) 


The approximations are valid as long as — A 7/2 ;§> 10 -I- 

rf). Alternatively we obtain for the case of an isothermal 
sphere with an on-axis shear (/? = 1, B = 2) 



2 2 
(1 -I- 7)2 (1 — 7)2(j2 


and 


ATr2 


-1 
(1 + 7 ) 



(1+7) \ 

(1 - iW ) 


( 20 ) 



(1 + 7)^ \ 
(l-7)Vy' 


• ( 21 ) 


Notice that the previous four equations are rotationally in¬ 
variant, and thus do not depend on whether we know the 
position angle of the lens galaxy. 


Based on numerical simulations, the two sets of equa¬ 
tions (eqs. 18 and 19) seem to be rather robust and they 
depend mainly on the uncertainty of D in eq. ( 8 ) which can 
be determined if the distances to the lensing galaxy and the 
source, i.e. the redshifts Zd, Zs and the Hubble constant Bq, 
are known. We note here that it is also possible to derive 
similar equations if we use the ratios of the time delay mo¬ 
ments or the image positions in the form of 


(1-7)" a %2 

EiW? ~ (1+ 7 )^'^ ~ Ar^2 


( 22 ) 


(see Appendix A). The equation is valid for the case 7 = 0 
or B = 2, /? = 1. Such equations have the advantage they are 
ratios and thus scale-free. However, unlike eqs. (18) to (21), 
to apply this relation, we need to transform the observed 
system into an unrotated coordinate system aligned with 
the symmetry axes of the lens. The reader should consult 
the Appendix C for the procedure to do this. 

Furthermore, when the external shear is negligible for 
an observed quadruple lens we can derive an approximate 
estimate on the slope p of the potential using the four image 
positions and the three time delays. The details can be found 
in the Appendix D. 


4 APPLICATION TO THE LENS SYSTEM 
RX J1131-1231 

In this section, we apply our formalism to the quadruple 
lens RX J1131—1231. The lensing galaxy is an elliptical at 
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Table 1. Image, galaxy positions and time delays for RX 
J1131—1231. The time delay tax is the time delay between image 
A and X = B,C, D. G denotes the lens galaxy. 


Object 

X 

[” 

1 

Y 

[” 

1 

tax [days] 

A 

-1-0.000 

± 

0.000 

+0.000 

± 

0.000 

0 

B 

-1-0.032 

± 

0.002 

+1.188 

± 

0.002 

0.7+ 1.0 

C 

-0.590 

± 

0.003 

-1.120 

± 

0.003 

0.0+ 1.3 

D 

-3.112 

± 

0.003 

+0.884 

± 

0.003 

93.8 + 1.5 

G 

-2.016 

± 

0.002 

+0.610 

± 

0.002 

- 


redshift 0.295 and the source redshift is 0.658. The large 
separation (« 3 arcsec) makes study of the system rela¬ 
tively simple (see Claeskens et al. 2006). The image, galaxy 
positions and the measured time delays (Tewes et al. 2013) 
are listed in Table 1 and Figure 1. For the cosmology we use 
a matter density of flm = 0.3, a cosmological constant of 
Ha = 0.7, and Hq = 70kms“^ Mpc“^. 

We first check whether the power-law potential gives a 
reasonable fit to the data. Witt (1996) showed that in such 
case, the images and the lensing galaxy are located on the 
same hyperbola-like curve. Furthermore, the position angle 
of the lens galaxy can also be inferred; the procedure to do 
this is outlined in the Appendix C. 

The results are shown in Fig. 1. Indeed the galaxy and 
images fall approximately on the curve, with the maximum 
distance from the hyperbola being about 0.025 arcsec. This 
should be compared with the positional accuracy of about a 
few mas (see Table 1). Our simple model appears to provide 
a reasonable approximation. Note that, from this model we 
can also derive the position angle of the lens galaxy, and 
obtain good agreement with the singular isothermal ellipsoid 
model in Claeskens et al. (2006). 

Recall that we implicitly assumed the angles are nor¬ 
malised by the Einstein radius, Oe (cf. eq. 4), i.e., a = 1. 
After normalisation, for /3 = 1,B = 2, we find for the sum 
of the squared image positions 



15.266(1")^ ^ 2 


(23) 


For the moment of the time delays in eq. (19) we need 
to relate the normalised units to the observed time delays 
and image positions. We find that D/(2c)(l -|- Zd) = 32.5 
days/(l”)^, and so 


-237.8 days (1")^ 
32.5 days 0'^ 



(24) 


From these two equations we find q = 0.69 and 9e = 1.57”. 
These values are roughly consistent with the results from the 
singular isothermal ellipsoid model in Table 4 in Claeskens 
et al. (2006), for which they gave 6e = 1.82”, and e = 0.45. 
The latter corresponds to an axis ratio q = 0.61 by the 
relationship 
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Figure 1. Image positions and galaxy position of RX 
J1131 —1231. The hyperbola line indicates the confinement of the 
galaxy and source positions if a pure elliptical potential (plus an 
on-axis shear) is assumed. The minimum distance of the galaxy 
position to the solid line is = 0.025” indicating that the 

lens system can not be perfectly fitted by any elliptical potential. 
The major axis (dashed line) is rotated by 6x/e — —74.63° which 
corresponds to a rotation angle of 6 = —15.37° from the r-axis 
(cf. Appendix C). 


a few percent enters from the uncertainty in the measured 
time delays. 


4.1 The impact of the slope of the potential on 
the model 

Such a simple elliptical isothermal model as discussed above 
can account only roughly the observed time delays. However, 
the model above would predict roughly AIa.s = —3 days, 
AIa.c = —19 days and AtA,D = 103 days using eq. (10) 
directly. This is certainly quite far away (~ lOcr) from the 
observed time delays as listed in Table 1. It even predicts 
that the order of arrival time would be C, B, A, D. 

In this section we will adjust the slope of the potential 
P in order to better account for the observed time delays, 
and see how the change affects the predicted axis ratio q. 
In the Appendix D we show how the slope of the potential 
can be inferred just from the observed image positions and 
time delays for a quadruple lens when the external shear 
is negligible. For the case of RX J1131—1231, we obtain 
B = 2.65 which corresponds to /3 = 1.14. 

Using eqs. (18) and (19) we have 



Our value differs from theirs only by approx. 10%. The dif¬ 
ferences are mostly due to the neglected source positions in 
eqs. (18) and (19) and the differences between elliptical den¬ 
sity and elliptical potential models. An additional error of 


1 + g^ 

1 - qB 


C 

y/B 


with 


C = 


1 

V-2AT,2 



3.99 (26) 


which depends entirely on the observed image positions and 
time delays of a quadruple lens. Finally we can deduce the 
axis ratio by 
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q 


c + Vb) 


0.72 


(27) 


This value differs by 5% of the previous derived value {q « 
0.69), and thus rather robust. For the Einstein radius we 
deduce 9e = 1.61” which differs about 3% from the previous 
model (9e = 1.57”). In comparison, this model with free 
slope predicts AIa.s ~ —1.1 days, AtA,c ~ —1.2 days and 
AtA,D ~ 93 days, closer to the observed values. Note that 
AtA,B and AtA,c depend rather sensitive on the goodness 
of the fit. Even this model requires further refinement, for 
example, the inclusion of an external shear will better fit the 
observed image positions. However, this may even change 
our derived axis ratio q by a few percent. 


4.2 Relations to the Hubble constant 


Finally we would like to point out a simple relation between 
the ellipticity e of the lensing galaxy and the Hubble con¬ 
stant Hq. Using eqs. (18), (19) and (25) we can write 

ATf 2 2c ATj ^2 

^ D{l + z,) 

A'T’ 

= ( 28 ) 

where Ri and ATe? denote the observed (unnormalised) 
quantities and 


f{Zd,Zs) 


2clHo 
71(1 -t- Zd) 


(29) 


is a function which just depends on the redshifts of the 
source and the lensing galaxy, the matter density and the 
cosmological constant, but not on the Hubble constant. 
However, this relation holds only if the lensing galaxy is 
not distorted by shear. 


5 SUMMARY AND DISCUSSION 

In the future, wide-held surveys such as the LSST will dis¬ 
cover tens of thousands of new gravitational lenses (Oguri & 
Marshall 2010), and for hundreds of these, their time delays 
will be measured accurately (Liao et al. 2015). Such sys¬ 
tems would be useful for independent determinations of the 
Hubble constant and understanding the evolution of dark 
energy. 

In this paper, we have developed a new approach by 
studying the cross-moments of the time delays and image 
positions with an elliptical power-law potential. We hnd that 
combinations of image positions and time delays for quadru¬ 
ple lenses yield simple analytical expressions that are con¬ 
nected with observable quantities, which can be used to ob¬ 
tain the approximate axis ratio q, the Einstein radius and 
the slope beta. We apply this method to RX J1131 —1231, 
and show that our analytical results match the full numeri¬ 
cal determinations approximately. 

These results may be useful for the following applica¬ 
tions. First, the derived axis ratio and slope can be used 
as initial guesses for further refinement through numerical 
modelling of the lenses. Second, the current Time Delay 
Challenge One (TDCl) (Liao et al. 2015) uses only the light 


curves, but not the image positions. In the next challenge, it 
is likely that the image positions will be taken into account 
in determining the time delay. Our moment approach may 
be useful in rejecting “catastrophic outliers” in the time de¬ 
lay. Third, our simple results may be useful for automated 
discovery of lenses in large surveys (Marshall et al. 2009; 
Chan et al. 2014). For example, the fact that the lens posi¬ 
tions fall on an approximate parabola may be useful in re¬ 
jecting unrealistic image configurations without doing much 
modelling. 

Our approach is clearly not accurate enough for each in¬ 
dividual system. However, if it is unbiased then it may still 
be useful for inferring the cosmological parameters (such 
as the Hubble constant Hq) statistically, similar to the ap¬ 
proach taken by (Oguri 2007). It will be interesting to study 
this issue using a mock sample of lenses in a future work. 
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APPENDIX A: MOMENTS OF THE TIME 
DELAY 

A1 The elliptical power-law potential 

For the elliptical potential given in eq. (1), the normalised 
time delay map is given by 

t{x, y) = {x- 0^ + {y-vf + y(Al) 

for 0 < /3 < 2. We show below that the time delay moments 
for /3 = 1 and /? = 4/3 can be derived exactly. These re¬ 
sults are then used to extrapolate to other values of /? with 
0 < /3 < 2. Note that we use a = 1 in this Appendix for 
convenience. 


A 1.1 Preliminaries 


Let us denote Ti = T{xi,yi) where {xi,yi) are the image 
positions of the n images i — if the the source is 

located inside the diamond caustic. Note that we have four 
images (n = 4) for 0 < ^ 1, while an additional central 

fifth image appears (n = 5) when 1 < /3 < 2. 

For the case /3 = 1 and /3 = 4/3 we can readily separate 
the solutions of the lens equation (x, y) into two polynomials. 


9^{x,i,ri) = ^aW , gy{y,i,'n) = ^biy , (A2) 

2=0 2=0 

with n = 4 for /? = 1 and n = 5 for j3 = 4/3 (cf. Witt & 
Mao 2000 and Hunter & Evans 2001). The moments of the 
image positions can now be related to the coefficients of the 
polynomials: 


= 


5,2 = 


S^3 = 


Q-n —1 

Q-n 


J2xi = - 
2 = 1 
n 


an-2 


(A3) 

(A4) 


yx? = -f^y + 3 ^"-"^"-^ -3^(A5) 

\ an / ai an 


o _ 4 (an-lY .^n-2an-l 


^^an — san—l ^an — 4 

an an an 


(A6) 


For the y component, analogous results are obtained: we just 
need to replace the ai by bi. For mixed terms we can use the 
relationship of the image positions 

- Xi)yi = Xiiri - yt) {A7) 


which is valid for any elliptical potential (cf. Witt 1996). 
Now we can write 


Sxy — ^ ^ Xiyi 
2 = 1 


C a _ 

(l_g2) y (l_g2)‘ 


Using these relations we are able to compute any desired 
mixed terms like Sxy 2 and so on successively. 

Recalling eq. (24) in Witt et al. (2000) , we can write 
for the time delay 

x{x, y) = i^ A + vy)- ^ {x^ + y^), (A9) 

where we expressed /? = 25/(2 -t- B) in terms of B. This 
equation is in particular valid because ip obeys the partial 
differential equation /3-i/) = xipx + y'4>v This enables us to 
express the potential tp in terms of the lens equation. 

For the moments of the time delay, we can write now 

n 

J2ri=n{e+V^)+^^^{^Sx+vSy)-^{Sx2+Sy2){AW) 
2 = 1 

and similarly, higher order moments can be written as 

yriXi = Sxif+ V^) + +ySxy) 

2 = 1 

-^{Sx3+Sxy2) (All) 

and so on. 


A1.2 Moment relations for elliptical power-law potentials 


In this section we present the results for the moments for 
the power-law potential. These results are only exact for the 
case /3 = 1 and /3 = 4/3. However, for different /3 values 
(0 < /3 < 2) numerical results differ only by a few percent 
or less if ^ and rj are small (the source is sufficiently well 
aligned with the line of sight). 

For convenience we define the following quantities 

and Q=iyi-qA (A12) 

Note that 5 = 2 for /3 = 1 and 5 = 4 for fS = 4/3. 

For the sum of the time delay we obtain in general 


-T- - ^ (-1 ^ ^ 2 . 2 2 

Q2^+Q2’l 

2=1 

and for the first order moments we end up with 
4 


Tx = 


5+4 + 


5 5Q2qS 


^ Q4 


Ty = -rj 


5 2 4q^ 

qB BqB + 5Q2 


+2^2 y2^_ sil 

^ ' g4 ' Q4 

For the second order moments we obtain 

4 2 55, . 

“ “5“^ ”^5"^^^^ “ BqBQ-i 


, r ^ 1 2 

”^^5^4 + 

2e"(e"-yY^2^2^2(i + y")^ 




g® 


(A13) 


(A14) 


(A15) 


(A8) 


(A16) 
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Ty2 


2q^ 


Bq^ 


(1 




_[4+2w 

I ^ o J gS 5Q4 


+ 


B ' 2 

2q\e-rff 

g® 


(A17) 


Now taking the results for the moments of the image posi¬ 
tions (which are only exact for the cases /3 = 0,1 and 4/3) 


Sx = 


(l+f)?"), 

(A18) 

Sy = 


■ (3 + 

(A19) 

Sx2 = 

B 0 
2 + (l + ^)C" 

.22 2 g^ 2 

(A20) 

Sy2 

2 , /I . 73, 2 2 g^ 2 . Sg"* 2 

gs + ( + 2 Q4 C + Q4 y , 

(A21) 

we obtain the relative 
quadrupole lenses 

moments for the time 

delay for 

AT. = 

Tx-To^, 


(A22) 

ATy = 

Ty-To^, 


(A23) 

A 7/2 = 

.-7- .-7- ^x2 

1x2 - '0^- 

_ 2 2 

B~^ BqB ’ 

(A24) 

ATy2 = 

r 

's/2 -10 — 

^ 2 2 

Bg 2 s + QqB ■ 

(A25) 


The latter two equations hold only if 77 ^ 1 are small. We 
did not expand the above expressions further because they 
follow a rather complicated pattern from which one does not 
gain much insight. However, it is interesting to note that 

AT,. = AT^2 + Ar ,2 « -1(^1 - (A26) 


and 

4 


S ^2 + Sy 2 — 2 -I- — -I- (1 -f + rf') 

q z 

■^g 2 ^ g 2 ^ qB (^27) 


for small ^ and rj. The expressions depend only on the dis¬ 
tances of the image positions to the galaxy centre and the 
relative time delays between the images, and so are rota- 
tionally invariant. In addition the ratios of the relative time 
delay moments satisfy 


A7),2 Sy2 


(A28) 


Such equations can be used to cross check our results. How¬ 
ever in this case we need to derive the coordinates of the 
unrotated system (see Appendix C). 


A 1.3 Moment relations for an isothermal sphere plus an 
on-axis shear 

For the case of an isothermal sphere plus an on axis shear 
(Witt & Mao 1997) we can derive the moments in a similar 
way as in Appendix A or by using the scaling properties 


of the lens equation in Appendix B. We just state here the 
approximate results for convenience. For the image position 
moments we have 


~ (1 -b 7)2 ~ (1 - 7 ) 2 g 2 

and for the time delay moments we obtain 

A.-,- -1 1 

A7;2 


and 


ATT ,2 


( 1 - 1 - 7 )^ (1-7)(1+ 7 )^'?^’ 


-1 


(l_7)3q4 


(A29) 


(A30) 


(A31) 


We note here that we only obtain rather simple results for 
the relative time delay moments for the case of an isothermal 
sphere when /3 = 1 and 7 = 0 . For completeness the results 
are as follows 


AT. = (A32) 

AT. = (A33) 

ATx2 = -1+4 - 10 c" 

g 2 

+Q-® [q^e - -q"v" + (A34) 

AT - zl±E_in!^ 

at;^ 10^2 

+Q"® -q* + 2q^rf - IQq'^frf] . (A35) 


APPENDIX B: THE INCLUSION OF AN 
ON-AXIS SHEAR 

We are able to transform the entire results we have obtained 
so far if the shear is aligned with the symmetry axis of the 
lens potential, i.e., an on-axis shear 7 . The lens equation in 
this case is given by 

^ = a; + 7 * — aa;(a:^ + (Blj 

V = y-qy —y)* +y q) (B2) 

(j 2 


We now show that these equations can be transformed into 
eqs. (2) and (3) by a simple transformation. If we denote the 
transformed quantities with a prime it yields 


C' 

^ W- ^ 

VI + 7 \/l - 7 

(B3) 

/ 

X 

= * 7/1 + 7 v = y\/i-7 

(B4) 

/ 

a 


(B5) 


whereas j3' = j3 remains unchanged. Applying the above 
transformation, we obtain 


= a;'-aV(*'" + y'Vg'")AA-i (B6) 

/ / O-'y' / a , 12 : /2n/3'/2-1 CTD-tA 

q = y —+y q T ^ , (B7) 

q/2 

identical to equations (2) and (3) if we drop the prime. We 
note that the magnification transforms like 

fi'i = Mi(l - 7^) 


(B 8 ) 
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and the time delay transforms like 


T = T — 




in 


(1 + 7 ) ( 1 - 7 ) 


(B9) 


where we have used eq. (9). Since we are only interested in 
the relative time delay between two images the additional 
terms in the transformation cancel out and we obtain 


Ar/j = Arij (BIO) 

This means that if we assume an elliptical power-law poten¬ 
tial the time delays between the images remain unchanged 
under the above transformation. Using only the time delays 
we would not be able to disentangle the contribution of the 
on-axis shear 7 and axis ratio q. 


APPENDIX C: THE ROTATION ANGLE OF 
THE LENSING GALAXY 

Often the positions of the images and the lensing galaxy are 
given relative to one image of the quadruple lens or the lens¬ 
ing galaxy. Also the usual RA and DEC coordinate systems 
may not be aligned with the symmetry axes of the lensing 
galaxy, and so that we can not apply the derived formulae 
in §3 directly. Since we assume an elliptical potential we are 
able to derive the rotation angle of the major axis due to 
eq. (7) of Witt (1996). 

Let Atti = ai be the right ascensions and ASi = di the 
declinations of the observed relative image positions relative 
to image 1 (image A) respectively. The rotation angle of the 
major axis is then given by 

(np\ — 2(02^2^34 -f 03^3^42 + 04^4^23) fT'I'l 

with kij = ttidj — ttjdi and i, j = 2,3,4 which corresponds to 
the images B, C and D. For RXJ 1131-1231 we obtain from 
Table 1 tan(20) = —0.5945 which correspond to a rotation 
angle of 6 = —15.37°. 

We note here that the rotation angle has an ambiguity 
of 90° since 6 is limited to —45° ^ 0 ^ 45°. However, this 
ambiguity can be easily resolved as long as the lensing galaxy 
can be observed (see Witt 1996 for further discussion). We 
finally arrive at the unrotated image positions by 

f Xi\ _ f cos 6 sin6» \ /fli - OG \ 

\yij ^ \-sm9 cose) \d,-dGj ^ 

where [acydo) is the relative galaxy position. Note that 9 
is measured relative to the x-axis. In observations the axis 
of the lensing galaxy is measured north-to-east ^n/e which 
is given by 0 n/e = —90° — 9 = —74.63° for this case. The 
rotation angle agrees quite well with the results of Claeskens 
et al. (2006) (Table 4) for their singular isothermal ellipsoid 
models with pure ellipticity and/or on-axis shear. 


APPENDIX D: ESTIMATION OF THE SLOPE 
OF THE LENS POTENTIAL 

When the external shear is negligible for an observed 
quadruple lens we can derive an approximate estimate on 
the slope p of the potential using just the 4 image posi¬ 
tions and the 3 time delays. Following eq. (25) in Witt et al. 
(2000), we can write for the time delay 


Ari.j = ^(r| - r|) -h - Xj) + nW - yj)], (Dl) 

where we expressed /3 in terms of B. Since the ratios of the 
time delays are independent of the distance parameters we 
can write Ati_ 2 /Ati _3 = Ari, 2 /Ari ,3 for example. Finally 
we use all the three time delay combinations to eliminate 
Tj or B. Now we obtain three equations in the form of 


ci^ -I- C 277 = 0 

{B — 2)0377 — 2ci = 0 
{B - 2)c 3^ -h 2 c 2 = 0 
where 


0 

Ati,2 

Ail ,3 

Ail ,4 

Xl 

X2 

X3 

Xi 

2 

2 

2 

2 

7^1 

7-2 

7-3 

ri 

1 

1 

1 

1 

0 

Ati,2 

Ail ,3 

Ail ,4 

yi 

7/2 

7/3 

7/4 

2 

2 

2 

2 

ri 

7-2 

7-3 

7-4 

1 

1 

1 

1 

0 

Ail ,2 

Ail ,3 

Ail ,4 

yi 

7/2 

7/3 

y4 

Xl 

X2 

X3 

Xi 

1 

1 

1 

1 


(D2) 


(03) 


are the determinant of different combinations of the time 
delays and the image positions. Further we recall that for a 
pure elliptical potential the relation 


inixiyj - Xjyi) + Xjyji^yi - Xiy) + Xiyi{xjri - iyj) = 0(D4) 


holds for i,j = 1, 2, 3, 4 and i j (Witt 1996). This enables 
us to eliminate ^ and y successively. Since we need here 
in principle only 2 image positions we may actually obtain 
for a complicated system different values of B from different 
combinations of image positions. We average over all 4 image 
positions to obtain an estimate of /3 


B = 2-h2—- 

616263 




63(6161 -I- 6262) 


where 






-1 

1 

-1 

1 


Xl 

372 

3:3 

Xi 

ei = 






3717/1 

3:27/2 

3:37/3 

3:47/4 


1 

1 

1 

1 


-1 

1 

-1 

1 


7/1 

7/2 

7/3 

y4 

62 = 






3717/1 

3727/2 

3:37/3 

3:47/4 


1 

1 

1 

1 


-1 

1 -1 

1 



7/1 

7/2 7/3 

7/4 


63 = 






371 

X2 X3 

Xi 



1 

1 1 

1 



(05) 


(06) 


The advantage of eq. (D5) is that it is completely scale free. 
We can immediately substitute the image positions and time 
delays without worrying about the units. However, we need 
to rotate the coordinate system by the angle of the lens 
galaxy 9 (see eq. Cl), which points to a disadvantage of 
eq. (D5) since the uncertainty on 9 enters the slope of the 
potential directly. 

For RX J1131 —1231 we obtain for a rotation angle of 
9g = —15.37°, B = 3.59 which correspond to /3 = 1.28. This 



9 


prediction seems to overestimate the slope of the potential 
because At ad is reduced to about At ad ~ 80 days. We 
can adjust the value of B so that the predicted time delays 
better match the observed values. We find that a value of 
B = 2.65 (/3 = 1.14) seems to be more appropriate which 
predicts At ad ~ 93 days comparable to the value in Table 
1. However, the large predicted offset is due to the neglect 
of an external shear. It would be desirable to extend eq.(D5) 
in such a way that it would include an off-axis shear. 

It is interesting that the system seems to be close to 
isothermal, but not exactly. This explains why we cannot 
account for the observed time delays very well. The conclu¬ 
sion is consistent with the more advanced model in Suyu 
et al. (2013) where they also had to allow a small deviation 
of the slope of the potential from isothermal (which corre¬ 
sponds to roughly p « 1.05 in our notation), but also include 
an external shear in their model. 
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ABSTRACT 

Time delays in gravitational lenses can be used to determine the Hubble constant and 
the lens potential. In future surveys, many gravitational lenses can be discovered, and 
their time delays and image positions can in principle be measured. Using an elliptical 
power-law potential, we show that combinations of image positions and time delays for 
quadruple lenses yield simple analytical expressions that are connected with observable 
quantities. These relations can be used to obtain the approximate axis ratio q, the 
Einstein radius and the slope. We apply this method to RX J1131—1231, and show 
that our analytical results match the full numerical determinations approximately. 
Our approach can quickly determine rough values of lens parameters, which can then 
be used as initial guesses for further refinement through numerical modelling and may 
be useful for automated lens search in large surveys. 

Key words: Gravitational leasing: strong - galaxies: structure - quasars: individual: 
RX J1131-1231 


1 INTRODUCTION 

Time delays in gravitational lenses provide an independent 
way of measuring the Hubble constant (?). Many studies 
have been performed in this direction. Observationally, two 
dozens or so lenses have their time delays accurately mea¬ 
sured after painstaking efforts (e.g., ?). Theoretical models 
also become increasingly sophisticated. For example, a re¬ 
cent study by ? takes into account of the line of sight struc¬ 
tures in addition to the lensing galaxy. Future surveys such 
as LSST will discover tens of thousands of lenses. In addi¬ 
tion, hundreds of lensed supernovae will be found. A large 
fraction (~ 40%) will have crude time delays (?), while a 
smaller fraction of these (~ 400) will have well-measured 
values (?). Such a large sample of lenses with time-delays 
can in principle be a very efficient way of constraining dark 
energy (?). 

The time delay for an isothermal potential or density 
distribution does not depend on the angular profile of the 
lensing galaxy (?). Other than this, very few analytical re¬ 
sults are known when the profile deviates from the isother¬ 
mal shape. With thousands of time delays becoming avail¬ 
able in the future, it is important to understand better how 
the time delays depend on the density prohle, and their im- 
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pact on the probability distribution of time delays. In this 
case one should not only focus on the time delays and the ra¬ 
tios of time delays measured in double and quadruple lenses. 
The ratio of time delays may depend on the lens model in a 
complex way. Therefore a diagram (of ratios) of time delays 
may be difficult to interpret. Many Monte Carlo simulations 
might be necessary to understand their behaviour and the 
dependency of different lens models on the parameters. In 
this paper we show that it can be useful to focus on the sum 
of time delays, and combinations of time delays and image 
positions. These quantities are more naturally connected to 
properties of the lensing galaxy such as the axis ratio or the 
slope of the inner potential. 


In this work we shall recall the properties of the ellipti¬ 
cal power law potential (§2). While this lens potential is not 
as realistic as the elliptical density distribution (e.g., ??), it 
often resembles the latter when the ellipticity is small (?) 
and provides useful insights into the lensing properties. In 
§3, we introduce a new moment approach through which we 
derive new connections between the time delays and their 
image positions. We apply the results to the quadruple lens 
RX J1131 —1231 in section §4, followed by a short discussion 
in §5. For ease of reading, we put most of the mathematical 
details in the appendices. 
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2 LENSING PROPERTIES OF AN 

ELLIPTICAL POWER-LAW POTENTIAL 

The elliptical power-law potential widely used in gravita¬ 
tional lensing (e.g., ??) is given by 

■fpix, y) = ^ for o < /3 < 2 (i) 

where {x, y) are the (angular) lens position, q is the axis ratio 
of the lensing galaxy (0 < g ^ 1), the slope of the potential 
and a sets the angular scale of the lens. For /3 = 1 we have the 
elliptical isothermal sphere which has very simple properties 
(see below). 

2.1 Lens equation and Einstein radius 

The lens equation for the power-law potential is 

C = x-^=x-ax{x^+y^ (2) 

^ ^ (3) 

where (^, r]) are the source position. If the lens is circularly 
symmetric (g = 1) and perfectly aligned with the source, the 
images form an Einstein ring with its angular scale given by 

= (4) 


2.2 Image Positions, Numbers and Magnifications 


The magnification of an image is generally given as 


1 _ f dr] dq\ ^ 

det J \dx dy dy dx J 


(5) 


where {xi,yi) is the image position of the ith image in the 
lens plane. The image positions where the magnifications 
are infinite form critical curves and their source positions 
form caustics. For an elliptical potential, the caustics are 
in general diamond-shaped. When a source falls inside the 
diamond caustics, we can have either four images for 0 < 
/3 ^ 1 or five images for 1 < /3 < 2. In the latter case, a 
faint, central image appears. 

The total magnification is given as the sum of all abso¬ 
lute magnifications 






i 


( 6 ) 


For an elliptical power-law potential the sum of the signed 
magnification follows a simple relation 




(7) 


for a quadruple system when the source is located inside the 
caustics. The result is exact for the cases /3 = 1 and p = 4/3 
(? and ?). 


2.3 Time Delay 

The time delay in gravitational lenses is given in general by 

(e.g. ?) 


( 8 ) 


T{x,y) = {x-pf+ {y-rif-2'4>{x,y), (9) 

where Zd is the redshift of the lens, and D = DdDs/Dds, 
with Dd,Ds and Dda are the angular-diameter distances to 
the lens, to the source, and from the lens to the source, 
respectively. Below we shall focus on r, and ignore all the 
cosmological dependences. 


3 TIME DELAYS AND IMAGE POSITIONS IN 

QUADRUPLE LENSES 

As mentioned before, when a source is located inside the 
caustics, we may observe at least 4 images and possibly 3 
time delays if the lensed source varies. Such observations of¬ 
fer an exciting opportunity to probe the (inner) potential of 
the lensing galaxy. In particular, future surveys will observe 
thousands of gravitational lensed systems (see ?). Some of 
theses systems might be fairly isolated with negligible shear, 
which, as we show below, would be an ideal case to study 
the potential slope and the axis ratio of the lensing galaxy 
in more details. In any case, the effects of a (small) shear 
enter linearly (see §A1.3), and so can be ignored as a first 
approximation. 

For an isothermal potential with arbitrary angular pro¬ 
files, 'i/(r,/)) = rF{(j)), the time delay between two images 
depend only on the distances of the images (xi, yi) from the 
galaxy centre, namely 

Ari,j = T{xi, yi) - T(xj,yj) = r| - r/ (10) 

with Vi = ytx'^ + yf (cf. ?). In particular for an isothermal 
elliptical potential it is rather straightforward to derive some 
invariants for the time delay. We can extend the results of 
? and obtain some simple relations concerning the moments 
of the magnification and image positions. Similar moment 
relations can be computed for the time delay as well. In 
particular one should not consider only the ratio of the time 
delays, but also the sum of all time delays of a system. We 
pursue this approach below; most of the mathematical de¬ 
tails are collected in the Appendices. 


3.1 A Moment Approach 


We will develop here a moment approach using the time de¬ 
lays to study the deviation from the ideal case of an elliptical 
isothermal sphere. Let us denote the moments of the time 
delay by 


= '^nxT, 

i 

7 - \ ^ m 

yrr. = Z^nyt , 


( 11 ) 

( 12 ) 


where n = rixi, yi) is the time delay at the i-th image posi¬ 
tion. Since the absolute time delays n can not be measured 
directly we relate them to measurable quantities by 


+n ^ xT 

(13) 

^Ari.iyr + Ti yT 

(14) 


i i 


Now we use the Oth-moment equation 


At 


— {l + Zd)T{x,y), 
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To = ^ Ari,! + nri (15) 


to eliminate the unknown time delay of the first image ri, 
where n is the total number of images, n is usually 4 (n = 5 
if a central image is present). Finally we can write 


AT.™ 

i 


and 

i i i 

(16) 

ATy^ 

= Tym - y 

n 

i 



= ^Ar,,iyr-i^Ar,,i^yr 

(17) 


i i i 


On the left side of the two previous equations we have now 
quantities which can be expressed for example in terms of 
rj, q and /3. On the right side, we have only terms involving 
the image positions and the relative time delays, which can 
be measured directly for a quadruple lens system. 


with smaller q. The more elliptical the system becomes the 
more the image positions are located further away from the 
Einstein ring. In particular in clusters where one expects 
1 < /3 < 1.5 this effect is considerably amplihed since for 
/3 = 4/3, 5 = 4 and for /? = 3/2, 5 = 6 . Therefore in 
clusters the source location plays less of a role and the sys¬ 
tem is mostly determined by q and j3. Similar effects can 
be observed if we have an isothermal sphere (/3 = 1 ) plus 
shear. The shear acts like an amplifier for smaller q. There¬ 
fore the shear 7 may emulate a more elliptical system than 
the ellipticity of the galaxy actually indicates. 

For the relative time delay moment A 7/2 the term be¬ 
comes small when the system is nearly round and q < 1. In 
this case the time delay is dominated by the source position 
and the approximation is fairly bad. For these cases we need 
to take the source positions into account. This would mean 
we need to solve a nonlinear set of equations and take the 
hrst order moments into account as well. Also in the case 
when the source is located near a cusp ^ or 7 might not 
be negligible any more. However, for rather elliptical sys¬ 
tem with q < 0.75 and /3 « 1 the source position plays less 
a role, so that the approximations hold fairly well in the 
case of galaxy leasing. Such a system would be suitable to 
determining the axis ratio q and the slope of the potential 


3.2 Analytical moment relations 


In the Appendix A we present the detailed procedure 
through which we find some analytic moment relations in¬ 
volving the relative time delay moments for a power-law po¬ 
tential. For simplicity, we normalise the coordinates by the 
Einstein radius, after which we have a = 1. For the second 
order moments we hnd, when the source is located close to 
the origin with no shear (7 = 0 ), 




with 


5 = 


2/3 

2-/3 


and 

A7;2 = IS.%2 +ATy2 


(18) 


(19) 


The approximations are valid as long as — A 7/2 10 (^^ -I- 

7 ^). Alternatively we obtain for the case of an isothermal 
sphere with an on-axis shear (/3 = 1, 5 = 2) 


2 ^ 2 

~ (1 + 7)" 

i=l 


(1 -7)2q2 


and 


AT /2 


-1 

(1 + 7) 



(1 + 7 ) \ 

(1 - 7 ) 52 ^ 


( 20 ) 



(1 + 7 )^ \ 
(l- 7 )Vy' 


• ( 21 ) 


Notice that the previous four equations are rotationally in¬ 
variant, and thus do not depend on whether we know the 
position angle of the lens galaxy. 

For fairly round system with g —>■ 1 the image positions 
are very close to the Einstein ring so that the solutions ap¬ 
proach ri —>■ 1 and the sum of the squared image positions 
^. r^ —>■ 4. The time delay between the images becomes 
smaller and smaller so that the relative time delay moment 
vanishes (AT /2 —>■ 0). For more elliptical systems the sum of 
the squared image distances to the galaxy centre increases 


P- 

Based on numerical simulations, the two sets of equa¬ 
tions (eqs. 18 and 19) seem to be rather robust and they 
depend mainly on the uncertainty of D in eq. ( 8 ) which can 
be determined if the distances to the lensing galaxy and the 
source, i.e. the redshifts Zd, Zs and the Hubble constant 5o, 
are known. We note here that it is also possible to derive 
similar equations if we use the ratios of the time delay mo¬ 
ments or the image positions in the form of 


B a%2 
EiW? ~ (1 + 7 )^'^ ~ ^Ty2 


( 22 ) 


(see Appendix A). The equation is valid for the case 7 = 0 
or 5 = 2, /3 = 1. Such equations have the advantage they are 
ratios and thus scale-free. However, unlike eqs. (18) to (21), 
to apply this relation, we need to transform the observed 
system into an unrotated coordinate system aligned with 
the symmetry axes of the lens. The reader should consult 
the Appendix C for the procedure to do this. 

Furthermore, when the external shear is negligible for 
an observed quadruple lens we can derive an approximate 
estimate on the slope /3 of the potential using the four image 
positions and the three time delays. The details can be found 
in the Appendix D. 


4 APPLICATION TO THE LENS SYSTEM 

RX J1131-1231 

In this section, we apply our formalism to the quadruple lens 
RX J1131—1231. The lensing galaxy is an elliptical at red- 
shift 0.295 and the source redshift is 0.658. The large separa¬ 
tion (« 3 arcsec) makes study of the system relatively simple 
(see ?). The image, galaxy positions and the measured time 
delays (?) are listed in Table 1 and Figure 1. For the cos¬ 
mology we use a matter density of Hm = 0.3, a cosmological 
constant of Ha = 0.7, and Ho = 70kms“^ Mpc“^. 

We first check whether the power-law potential gives 
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Table 1. Image, galaxy positions and time delays for RX 
J1131—1231. The time delay tax is the time delay between image 
A and X = B,C, D. G denotes the lens galaxy. 


Object 

X 

[” 

1 

Y 

[” 

1 

tax [days] 

A 

-1-0.000 

± 

0.000 

+0.000 

± 

0.000 

0 

B 

-1-0.032 

± 

0.002 

+1.188 

± 

0.002 

0.7+ 1.0 

C 

-0.590 

± 

0.003 

-1.120 

± 

0.003 

0.0+ 1.3 

D 

-3.112 

± 

0.003 

+0.884 

± 

0.003 

93.8 + 1.5 

G 

-2.016 

± 

0.002 

+0.610 

± 

0.002 

- 


a reasonable fit to the data. ? showed that in such case, 
the images and the lensing galaxy are located on the same 
hyperbola-like curve. Furthermore, the position angle of the 
lens galaxy can also be inferred; the procedure to do this is 
outlined in the Appendix C. 

The results are shown in Fig. 1. Indeed the galaxy and 
images fall approximately on the curve, with the maximum 
distance from the hyperbola being about 0.025 arcsec. This 
should be compared with the positional accuracy of about a 
few mas (see Table 1). Our simple model appears to provide 
a reasonable approximation. Note that, from this model we 
can also derive the position angle of the lens galaxy, and 
obtain good agreement with the singular isothermal ellipsoid 
model in ?. 

Recall that we implicitly assumed the angles are nor¬ 
malised by the Einstein radius, Oe (cf. eq. 4), i.e., a — 1. 
After normalisation, for /3 = 1,B = 2, we find for the sum 
of the squared image positions 



15.266(1")^ ^ 2 

el 


(23) 


For the moment of the time delays in eq. (19) we need 
to relate the normalised units to the observed time delays 
and image positions. We find that D/(2c)(l -|- Zd) = 32.5 
days/(l")^, and so 


-237.8 days (1")^ 
32.5 days 9^^ 



(24) 


From these two equations we find q = 0.69 and 9e = 1.57”. 
These values are roughly consistent with the results from the 
singular isothermal ellipsoid model in Table 4 in ?, for which 
they gave 9e = 1.82”, and e = 0.45. The latter corresponds 
to an axis ratio q = 0.61 by the relationship 



Our value differs from theirs only by approx. 10%. The dif¬ 
ferences are mostly due to the neglected source positions in 
eqs. (18) and (19) and the differences between elliptical den¬ 
sity and elliptical potential models. An additional error of 
a few percent enters from the uncertainty in the measured 
time delays. 


4.1 The impact of the slope of the potential on 
the model 

Such a simple elliptical isothermal model as discussed above 
can account only roughly the observed time delays. However, 



Aa[arcseo] 

Figure 1. Image positions and galaxy position of RX 
J1131 —1231. The hyperbola line indicates the confinement of the 
galaxy and source positions if a pure elliptical potential (plus an 
on-axis shear) is assumed. The minimum distance of the galaxy 
position to the solid line is = 0.025” indicating that the 

lens system can not be perfectly fitted by any elliptical potential. 
The major axis (dashed line) is rotated by 9x/e — —74.63° which 
corresponds to a rotation angle of 6 = —15.37° from the r-axis 
(cf. Appendix C). 


the model above would predict roughly AIa.s = —3 days, 
AtA.c = —19 days and AtA,D = 103 days using eq. (10) 
directly. This is certainly quite far away (~ lOcr) from the 
observed time delays as listed in Table 1. It even predicts 
that the order of arrival time would be C, B, A, D. 

In this section we will adjust the slope of the potential 
P in order to better account for the observed time delays, 
and see how the change affects the predicted axis ratio q. 
In the Appendix D we show how the slope of the potential 
can be inferred just from the observed image positions and 
time delays for a quadruple lens when the external shear 
is negligible. For the case of RX J1131—1231, we obtain 
B = 2.65 which corresponds to /3 = 1.14. 

Using eqs. (18) and (19) we have 


1-bg® 
1 - qB 


C 

y/B 


with 


C = 


1 

V-2AT,2 



3.99 (26) 


which depends entirely on the observed image positions and 
time delays of a quadruple lens. Finally we can deduce the 
axis ratio by 


<1 


( C-^ \ 

\c + ^/b) 


l/B 


0.72 


(27) 


This value differs by 5% of the previous derived value [q « 
0.69), and thus rather robust. For the Einstein radius we 
deduce 9e = 1.61” which differs about 3% from the previous 
model (Oe = 1.57”). In comparison, this model with free 
slope predicts AIa.b ~ —1.1 days, AIa.c ~ —1.2 days and 
AtA,D ~ 93 days, closer to the observed values. Note that 
AtA,B and AtA,c depend rather sensitive on the goodness 
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of the fit. Even this model requires further refinement, for 
example, the inclusion of an external shear will better fit the 
observed image positions. However, this may even change 
our derived axis ratio g by a few percent. 


be useful for inferring the cosmological parameters (such 
as the Hubble constant Hq) statistically, similar to the ap¬ 
proach taken by (?). It will be interesting to study this issue 
using a mock sample of lenses in a future work. 


4.2 Relations to the Hubble constant 


Finally we would like to point out a simple relation between 
the ellipticity e of the lensing galaxy and the Hubble con¬ 
stant Hq. Using eqs. (18), (19) and (25) we can write 


A7).2 


2c ATji2 


(E.-'i: 


AT„ 


(E.BO 


rHof{zd,Zs 


(28) 


where Ri and ATji 2 denote the observed (unnormalised) 
quantities and 


f{Zd,Zs) 


2c/Ho 
D{1 Zd) 


(29) 


is a function which just depends on the redshifts of the 
source and the lensing galaxy, the matter density and the 
cosmological constant, but not on the Hubble constant. 
However, this relation holds only if the lensing galaxy is 
not distorted by shear. 


5 SUMMARY AND DISCUSSION 

In the future, wide-field surveys such as the LSST will dis¬ 
cover tens of thousands of new gravitational lenses (?), and 
for hundreds of these, their time delays will be measured ac¬ 
curately (?). Such systems would be useful for independent 
determinations of the Hubble constant and understanding 
the evolution of dark energy. 

In this paper, we have developed a new approach by 
studying the cross-moments of the time delays and image 
positions with an elliptical power-law potential. We find that 
combinations of image positions and time delays for quadru¬ 
ple lenses yield simple analytical expressions that are con¬ 
nected with observable quantities, which can be used to ob¬ 
tain the approximate axis ratio q, the Einstein radius and 
the slope beta. We apply this method to RX J1131 —1231, 
and show that our analytical results match the full numeri¬ 
cal determinations approximately. 

These results may be useful for the following applica¬ 
tions. First, the derived axis ratio and slope can be used 
as initial guesses for further refinement through numerical 
modelling of the lenses. Second, the current Time Delay 
Challenge One (TDCl) (?) uses only the light curves, but 
not the image positions. In the next challenge, it is likely 
that the image positions will be taken into account in deter¬ 
mining the time delay. Our moment approach may be useful 
in rejecting “catastrophic outliers” in the time delay. Third, 
our simple results may be useful for automated discovery 
of lenses in large surveys (?; ?). For example, the fact that 
the lens positions fall on an approximate parabola may be 
useful in rejecting unrealistic image configurations without 
doing much modelling. 

Our approach is clearly not accurate enough for each in¬ 
dividual system. However, if it is unbiased then it may still 
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APPENDIX A: MOMENTS OF THE TIME 
DELAY 

A1 The elliptical power-law potential 

For the elliptical potential given in eq. (1), the normalised 
time delay map is given by 

T{x,y) = {x-0^ -b {y-nf - ‘^{x^ +y^I (Al) 

for 0 < /? < 2. We show below that the time delay moments 
for /3 = 1 and p = A/?) can be derived exactly. These re¬ 
sults are then used to extrapolate to other values of /3 with 
0 < /3 < 2. Note that we use a = 1 in this Appendix for 
convenience. 


A 1.1 Preliminaries 

Let us denote Ti = T{xi,yi) where {xi,yi) are the image 
positions of the n images i = l,...,n if the the source is 
located inside the diamond caustic. Note that we have four 
images (n = 4) for 0 < /3 ^ 1, while an additional central 
fifth image appears (n = 5) when 1 < /3 < 2. 

For the case /3 = 1 and p = 4/3 we can readily separate 
the solutions of the lens equation {x, y) into two polynomials. 

n n 

g^ix,p,v) , gy{y,^,v) = '^b^y'^ , (A2) 

4 = 0 4 = 0 

with n = 4 for p = 1 and n = 5 for P = 4/3 (cf. ? and ?). 
The moments of the image positions can now be related to 
the coefficients of the polynomials: 

& SE 

S^2 := 

•5.4 = 


For the y component, analogous results are obtained: we just 


J2x^ = - 


O^n—l 

Q-n 


^Cln-2 

(^n 


(A3) 

(A4) 


4=1 

n 2 

y4 = -f —1 _3gEzi(A5) 

^ V a„ / al a„ ^ 

4=1 

n 42 

E 4 _ / ^n—1 \ ^ —2^n—1 

(A6) 
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need to replace the ai by bi . For mixed terms we can use the 
relationship of the image positions 

- Xi)yi = Xiiri - yi) (A7) 

which is valid for any elliptical potential (cf. ?). Now we can 
write 

n 2 

Sxy = ^ Xiyi = ~ (1 _ 

i=l 

Using these relations we are able to compute any desired 
mixed terms like S^y 2 and so on successively. 

Recalling eq. (24) in ? , we can write for the time delay 

Tix,y) =f+y'^+ {^x + yy) - ^{x^ +y'^), (A9) 

where we expressed /3 = 25/(2 + B) in terms of B. This 
equation is in particular valid because ip obeys the partial 
differential equation piip = xip^ + yipy This enables us to 
express the potential ip in terms of the lens equation. 

For the moments of the time delay, we can write now 

{^S^+ySy)-^{S^2+Sy2){A10) 

i=l 

and similarly, higher order moments can be written as 

n . . 

'^TiXi = + 77 ^) + — {^S^2+yS^y) 

i=l 

-^{S^3+S^y2) (All) 

and so on. 


A 1.2 Moment relations for elliptical power-law potentials 


In this section we present the results for the moments for 
the power-law potential. These results are only exact for the 
case /? = 1 and jS = 4/3. However, for different /3 values 
(0 < /3 < 2) numerical results differ only by a few percent 
or less if ^ and y are small (the source is sufficiently well 
aligned with the line of sight). 

For convenience we define the following quantities 

B=i^ and Q=^l-q2. (A12) 

Note that B — 2 for /3 = 1 and 5 = 4 for /3 = 4/3. 

For the sum of the time delay we obtain in general 


ro = ^T, 

i = l 




2^-2 . ‘2 2 


(A13) 


and for the first order moments we end up with 


r. = 


Ty 


2 4 

5 ^ 5Q2gS 


0^-3 , Of-,^2 (1 + ) 

- 2 ^ ^ + 2^77 


-V 


2 ^ 4g^ 


qB ^ BqB ' BQ2 




(A14) 


(A15) 


For the second order moments we obtain 

4C" 


7/2 = 


-[4-1-1-[C"’ _ 

5 *52^^ BqBQi 


+ [ 


+ 

2q\e-riW 




^2^2^2(i + g")" 


g® 


(AI 6 ) 


7/2 = — 


2 ?" 


Bq'^B (l-q2)2 


r 


2 55, 77^ Aq^y^ 


B^ 2\b bQ* 

2 q\e-vW . 2 . 22(1 + ?")" 

+-Qs- 2qiy^—. (A17) 

Now taking the results for the moments of the image posi¬ 
tions (which are only exact for the cases /3 = 0,1 and 4/3) 

= 

Sy = 

S ^2 

Sy2 = 

we obtain the relative moments for the time delay for 
quadrupole lenses 


^((3+|)-(l + |)?"), 

(A18) 

^((l + f)-(3+f)?"), 

(A19) 

9 , /I , -B, 2 , 2 2 2(7^ 2 

2 + (i+ 3 )? -^v, 

(A20) 

2 , , B. 2 2q^ 2 . 2 g'* 2 

(A21) 


A% 

11 

1 

5) 


(A22) 

ATy 

II 

1 

5) 


(A23) 

A%2 

9 - T B^2 

= /a:2 - '0^- - 

2 2 

5 BqB ’ 

(A24) 

ATy2 

= 5,2-To—S 

2 2 

Bq^B + BqB ’ 

(A25) 


The latter two equations hold only if 77 ^ 1 are small. We 
did not expand the above expressions further because they 
follow a rather complicated pattern from which one does not 
gain much insight. However, it is interesting to note that 

Ar.2 = Ar.2 + Ar,2«-|(^i-^^ (A26) 


and 

4 


S^2 -\- Sy2 —2-1-^ -I- (1 + — )(C" + ??") 

q z 

, _^.2 _ 2 9 , 7A97'\ 

■^g2^ Q 2 7 + qB ("^27) 


for small ^ and y. The expressions depend only on the dis¬ 
tances of the image positions to the galaxy centre and the 
relative time delays between the images, and so are rota- 
tionally invariant. In addition the ratios of the relative time 
delay moments satisfy 


A7/,2 ^ B ^ S^2 


(A28) 


Such equations can be used to cross check our results. How¬ 
ever in this case we need to derive the coordinates of the 
unrotated system (see Appendix C). 
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Al.S Moment relations for an isothermal sphere plus an 
on-axis shear 

For the case of an isothermal sphere plus an on axis shear (?) 
we can derive the moments in a similar way as in Appendix 
A or by using the scaling properties of the lens equation in 
Appendix B. We just state here the approximate results for 
convenience. For the image position moments we have 


~ (1 + 7)2 ~ (1 - 7 ) 2 q 2 

and for the time delay moments we obtain 


AX. 


and 


AX. 


-1 


+ 


1 


(1 + 7)3 (l-7)(l + 7)2g2’ 


-1 


1 


(l_7)3q,4 (1 _ 


(A29) 


(A30) 


(A31) 


We note here that we only obtain rather simple results for 
the relative time delay moments for the case of an isothermal 
sphere when /? = 1 and 7 = 0. For completeness the results 
are as follows 


AX = 

ax 2 = -1 + 4-ioe" 


(A32) 

(A33) 


+Q-® -q^n+ (A34) 


AT - 

qi 10 g2 


+Q-® [q*e -n+ - lOq^ev"] ■ (A35) 


APPENDIX B: THE INCLUSION OF AN 
ON-AXIS SHEAR 

We are able to transform the entire results we have obtained 
so far if the shear is aligned with the symmetry axis of the 
lens potential, i.e., an on-axis shear 7 . The lens equation in 
this case is given by 


^ = a: + 7 ® — ax{x'^ + ^ 

/ 2 , 2 / 2\p/2-l 

1 = y-iy - + 2 / /<? ) 


(Bl) 

(B2) 


We now show that these equations can be transformed into 
eqs. (2) and (3) by a simple transformation. If we denote the 
transformed quantities with a prime it yields 


^ \/l + 7 ^ VI -7 

x' = * 1/1 + 7 y' = 2 /\/i -7 

a' = a(l + 7)-'3/2 = 


(B3) 

(B4) 

(B5) 


whereas /3' = p remains unchanged. Applying the above 
transformation, we obtain 


/ / // /2 I /2 / i2\8'f2 — l 

4 = X — a X [x + y Iq ) 

2 V 

„I2 




(B 6 ) 

(B7) 


identical to equations (2) and (3) if we drop the prime. We 
note that the magnification transforms like 


y'i = - 7 / 

and the time delay transforms like 


+2 2 

74 7 ^ 

■ X — 7 :^-— 7 + 


( 1 + 7 ) ( 1 - 7 ) 


(B 8 ) 


(B9) 


where we have used eq. (9). Since we are only interested in 
the relative time delay between two images the additional 
terms in the transformation cancel out and we obtain 


At/ = Ari,. 


(BIO) 


This means that if we assume an elliptical power-law poten¬ 
tial the time delays between the images remain unchanged 
under the above transformation. Using only the time delays 
we would not be able to disentangle the contribution of the 
on-axis shear 7 and axis ratio q. 


APPENDIX C: THE ROTATION ANGLE OF 
THE LENSING GALAXY 

Often the positions of the images and the lensing galaxy are 
given relative to one image of the quadruple lens or the lens¬ 
ing galaxy. Also the usual RA and DEC coordinate systems 
may not be aligned with the symmetry axes of the lensing 
galaxy, and so that we can not apply the derived formulae 
in §3 directly. Since we assume an elliptical potential we are 
able to derive the rotation angle of the major axis due to 
eq. (7) of ?. 

Let Aa; = ai be the right ascensions and A(5i = di the 
declinations of the observed relative image positions relative 
to image 1 (image A) respectively. The rotation angle of the 
major axis is then given by 


tan( 20 ) = 


2(02X^34 + 03X^42 + a 4 d 4 , k 23 ) 

{al — d\)k34 + (o| — d|)fc42 + (a| — d|)fc23 


(Cl) 


with kij = aidj —ajdi and i,j = 2, 3, 4 which corresponds to 
the images B, C and D. For RXJ 1131-1231 we obtain from 
Table 1 tan(20) = —0.5945 which correspond to a rotation 
angle oi 9 = —15.37°. 

We note here that the rotation angle has an ambiguity 
of 90° since 9 is limited to —45° ^ 0 ^ 45°. However, this 
ambiguity can be easily resolved as long as the lensing galaxy 
can be observed (see ? for further discussion). We finally 
arrive at the unrotated image positions by 


cos 9 sin 9 
— sin 9 cos 9 


ai - aa 
di - da 


(C2) 


where {aajda) is the relative galaxy position. Note that 9 
is measured relative to the x-axis. In observations the axis 
of the lensing galaxy is measured north-to-east 0 n/e which 
is given by 0 n/b = —90° — 9 — —74.63° for this case. The 
rotation angle agrees quite well with the results of ? (Table 
4) for their singular isothermal ellipsoid models with pure 
ellipticity and/or on-axis shear. 


1 
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APPENDIX D: ESTIMATION OF THE SLOPE 
OF THE LENS POTENTIAL 


When the external shear is negligible for an observed 
quadruple lens we can derive an approximate estimate on 
the slope /? of the potential using just the 4 image positions 
and the 3 time delays. Following eq. (25) in ?, we can write 
for the time delay 

- Vi) + + VW - Vi)]’ 

where we expressed /3 in terms of B. Since the ratios of the 
time delays are independent of the distance parameters we 
can write Ati, 2 /Ati ,3 = Ari, 2 /Ari _3 for example. Finally 
we use all the three time delay combinations to eliminate 
r] or B. Now we obtain three equations in the form of 


ClC + 0211 = 

0 


{B - 2)csr) - 

2 ci = 

0 


[B - 2)c3e + 2 c2 = 

0 


where 

0 

Atx ,2 

Atl,3 

Atx, 


Xx 

X2 

X3 

X4 

61 = 

2 

2 

2 

2 


61 

62 

xs 

Xi 


1 

1 

1 

1 


0 

Atx ,2 

Atl,3 

Atx, 

62 = 

yi 

92 

93 

94, 

2 

2 

2 

2 



62 

r3 

Xi 


1 

1 

1 

1 


0 

Atx ,2 

Atx,3 

Atx, 

63 = 

Vi 

V2 

93 

94 


Xx 

X2 

X3 

X4 


1 

1 

1 

1 


(D2) 


(D3) 


are the determinant of different combinations of the time 
delays and the image positions. Further we recall that for a 
pure elliptical potential the relation 


^vixiVj - XjVi) + XjVji^yi - XiTi) + Xiyiixj-n - ^y/) = 0(D4) 


holds for i,j = 1,2, 3, 4 and i ^ j (?). This enables us to 
eliminate ^ and rj successively. Since we need here in prin¬ 
ciple only 2 image positions we may actually obtain for a 
complicated system different values of B from different com¬ 
binations of image positions. We average over all 4 image 
positions to obtain an estimate of /3 


B^2 + 2 
where 

Cl = 

62 = 


Cl 02 63 

C3(ciei - I - C 262 ) 


-1 

1 

-1 

1 

Xx 

X2 

X 3 

X 4 

X191 

X292 

X393 

X4y4 

1 

1 

1 

1 

-1 

1 

-1 

1 

91 

92 

93 

94 

X191 

X292 

X393 

X4y4 

1 

1 

1 

1 

-1 

1 

-1 

1 


94 

92 

93 

94 


Xx 

X2 

X3 

X4 


1 

1 

1 

1 



(D5) 


The advantage of eq. (D5) is that it is completely scale free. 
We can immediately substitute the image positions and time 
delays without worrying about the units. However, we need 
to rotate the coordinate system by the angle of the lens 
galaxy 6 (see eq. Cl), which points to a disadvantage of 
eq. (D5) since the uncertainty on 9 enters the slope of the 
potential directly. 

For RX J1131 —1231 we obtain for a rotation angle of 
9g = —15.37°, B = 3.59 which correspond to /3 = 1.28. This 
prediction seems to overestimate the slope of the potential 
because At ad is reduced to about At ad ~ 80 days. We 
can adjust the value of B so that the predicted time delays 
better match the observed values. We find that a value of 
B = 2.65 (/3 = 1.14) seems to be more appropriate which 
predicts At ad ~ 93 days comparable to the value in Table 
1. However, the large predicted offset is due to the neglect 
of an external shear. It would be desirable to extend eq.(D5) 
in such a way that it would include an off-axis shear. 

It is interesting that the system seems to be close to 
isothermal, but not exactly. This explains why we cannot ac¬ 
count for the observed time delays very well. The conclusion 
is consistent with the more advanced model in ? where they 
also had to allow a small deviation of the slope of the poten¬ 
tial from isothermal (which corresponds to roughly /3 « 1.05 
in our notation), but also include an external shear in their 
model. 


63 


(D 6 ) 



